An Eccentric Binary Millisecond Pulsar in the Galactic Plane 
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Binary pulsar systems are superb probes of stellar and binary evolution and 
the physics of extreme environments. In a survey with the Arecibo telescope, 
we have found PSR J1903+0327, a radio pulsar with a rotational period of 
2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass com- 
panion. Infrared observations identify a possible main-sequence companion 
star. Conventional binary stellar evolution models predict neither large orbital 
eccentricities nor main-sequence companions around millisecond pulsars. Al- 
ternative formation scenarios involve recycling a neutron star in a globular 
cluster then ejecting it into the Galactic disk or membership in a hierarchical 
triple system. A relativistic analysis of timing observations of the pulsar finds 
its mass to be 1.74±O.O4M , an unusually high value. 

The population of binary and millisecond pulsars in the disk of our Galaxy is thought to 
have two main formation mechanisms CO). Most pulsars with spin periods of tens of millisec- 
onds have neutron star (NS) companions in orbits of high eccentricity caused by the near dis- 
ruption of the system by a supernova explosion. In contrast, pulsars with spin periods less than 
about 10 ms (i.e. "millisecond pulsars" or MSPs) have white dwarf (WD) companions in orbits 
made highly circular (orbital eccentricities e < 0.001) by tidal effects during the recycling pro- 
cess. The combination of rapid spin rates and circular orbits is considered vital evidence that 
MSPs achieve their short periods via accretion of mass and angular momentum from binary 
companion stars ©. Here we report the discovery of an unprecedented MSP that requires a 
different formation mechanism and whose potentially large mass may play an important role in 
constraining the equation of state of matter at supra-nuclear density ©. 
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Discovery and Follow-up Observations 

We are conducting a pulsar survey of the Galactic plane using the Arecibo L-band Feed Array 
(ALFA) receiver on the 305-m Arecibo radio telescope in Puerto Rico (4). The large collecting 
area of Arecibo, the rapid sampling rate (every 64 /is), and the high spectral resolution (256 
channels over 100 MHz, which minimizes the dispersive smearing due to free electrons along 
the line of sight) provides sensitivity to MSPs over a much larger volume of the Galactic disk 
than any previous pulsar survey. 

The 2. 15-ms pulsar J 1903+0327 was discovered using a search pipeline based on the PRESTO 
suite of pulsar analysis software £5] It was detected as a highly significant signal with a large 
dispersion measure (DM) of 297pccm~ 3 in data taken in October 2005. Follow-up timing 
observations using Arecibo, the Green Bank Telescope, and the Westerbork Synthesis Radio 
Telescope revealed the binary orbit of the pulsar to be highly eccentric (Fig. 1). The Keplerian 
orbital parameters give a minimum companion mass of 0.85 — 1.07 M (for pulsar masses of 
1.3 — 1.9 M Q ). Additional constraints come from the extensive timing of the pulsar, including 
a general relativistic interpretation of the measurement of the advance of periastron, to, and a 
detection of the Shapiro delay. As described in the supplementary online material, and in the 
caption to Table 1 , the best model fit to the pulse times of arrival indicates that the companion 
has a mass of 1.051(15) M and that the pulsar has a mass of 1.74(4) M Q . It should be noted 
that these masses are based on ~1.5 years of timing data which cannot yet include parame- 
ters such as proper motion which can affect the mass measurements. While the pulsar mass is 
significantly more massive than the 1.25 — 1.45 M Q seen in most double neutron star (DNS) sys- 
tems © it is comparable to the inferred masses of several recently detected pulsars in eccentric 
binaries in globular clusters (0 |U |9|), at least one other Galactic MSP (iTJ), and the X-ray pulsar 
Vela X-l (|771) . If the large pulsar mass is confirmed in future observations, it will constrain the 
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equation of state of matter at supra-nuclear density and potentially rule-out certain 'soft' equa- 
tions of state (3). The companion mass is compatible with those of a NS, WD, or main-sequence 
(MS) companion. Although the spin parameters of PSR J1903+0327 (Table 1) resemble those 
of other Galactic disk MSPs, the pulsar is clearly distinct when orbital eccentricity is compared 
as well (Fig. 2). Further details of the search pipeline, follow-up observations, and the timing 
analysis are provided in the supporting online material. 

Given the possibility that the companion could be a NS, and potentially a pulsar, we searched 
for pulsations from the companion using several of the Arecibo observations, but found none. 
These null results set an upper limit on pulsed emission at 1.4 GHz of ~20 jiiy for a period of 
2 ms and ~9 fiJy for a period of 200 ms assuming a pulse duty cycle of 30% of the pulse period. 

To search for a MS companion, we obtained images of the pulsar field with the Gemini North 
telescope on July 24, 2007. The total exposure times were lOmin in the infrared J, H and K s 
bands (1.27, 1.67 and 2.22 /zm, respectively). After calibrating the astrometry and photometry 
of the images against the 2MASS catalog, we find a single star within the 0.13"1-ct frame-tie 
error circle at the position of the pulsar (Fig. 3). It has J = 19.22(9), H = 18.41(10) and 
Ks = 18.03(9) magnitudes. Given the density of stars in this field, we estimate the probability 
of finding a star in the error circle by chance is 2.6%. Using MS star models (1721) and estimating 
the reddening with red clump stars ((73)) at the ~6.4kpc distance inferred from the pulsar's DM 
(|74)) . we find that a 0.9 M star of age 10 Gyr would have similar magnitudes to those observed. 
The uncertainties in the distance and reddening measurements also allow for a 1.05 M star of 
age 1-5 Gyr making it possible that this (likely) MS star is a companion to PSR J1903+0327. 

If the MS star is the pulsar's binary companion, ionized stellar winds might be detected as 
increases in the measured DM or as an eclipse of the pulsed flux when the companion passes 
between the pulsar and Earth. From multi-frequency observations around the orbit, we find 
no evidence of an eclipse, and we limit any additional DM contribution to <0.02pccm~ 3 . If 
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the companion of PSR J 1903+0327 is a ~1 M main-sequence star with Solar- like winds, 
we would expect an additional DM contribution of order 10~ 3 pccm 3 near conjunction (1751) . 
Strong irradiation of the companion by the pulsar's relativistic wind, however, could lead to a 
substantially larger mass loss. Our DM variation limit argues against such a large mass loss. 

For a system like PSR J1903+0327, a rotationally induced quadrupole moment in the com- 
panion could cause a classical periastron advance (|7oT) that would contribute to the measured 
to = 2.46(2) x 10~ 4 degyr _1 . If the companion star is a WD rotating near breakup velocity, 
the classical contribution would typically be of order 10~ 7 degyr _1 , but could be more than 
an order of magnitude larger for specific but unlikely system orientations. If the companion 
is a 1— 2Gyr, 1.05 M Q main-sequence star with a typical rotational period of 8—10 days (1771) . 
the classical contribution would be 2x 10~ 6 degyr -1 for most system orientations. Such a star 
would need a rotational period between 1.3 and 1.5 days and/or an unlikely system orientation 
to account for ~10% of the measured to, or a rotational period <0.5 day to account for all it. 
These numbers suggest that the measured cu is dominated by general relativistic effects and, 
given the high quality of the fit, that the use of the relativistic timing model is well justified. 

Formation Mechanisms 

What is the origin of this unique system with a short spin period, large orbital eccentricity, and 
possible MS companion? According to conventional evolutionary scenarios CO), binary pulsars 
that have been recycled down to millisecond periods should always appear in circularized orbits. 
In contrast, pulsars in eccentric systems should be only mildly recycled or not recycled at all. 
Since PSR J1903+0327 does not fall into either of these broad categories, we consider three 
alternative scenarios. The first is that the pulsar was not recycled but was born spinning rapidly 
in an eccentric orbit at the time the NS was created. The second has the pulsar recycled in a 
globular cluster and then ejected into the Galactic disk. The third has the pulsar recycled in a 
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hierarchical triple system. 

For various reasons it seems unlikely that the pulsar was formed spinning rapidly at the time 
of core-collapse with a small surface magnetic field (2x 10 8 G). First, there are no pulsars like 
J1903+0327 in any of the more than 50 young supernova remnants in which a NS has been 
inferred or detected directly ( Tffi) . Second, a "born-fast" scenario for PSR J 1903+0327 would 
likely be expected to account for some or all of the 18 isolated MSPs detected in the Galactic 
disk. While the formation of those systems is puzzling since the observed timescales for evapo- 
rating a companion star seem too long, their spin distributions, space velocities, and energetics 
are indistinguishable from those of recycled (i.e. not born-fast) binary MSPs and their space 
velocities and scale heights do not match those of non-recycled pulsars (17"9l) . Third, magnetic 
fields in young pulsars likely originate either from dynamo action in the proto-NS (|20l) or via 
compression of 'frozen-in' fields of the progenitor star during collapse (1271) . If compression is 
the correct mechanism, then young pulsars with magnetic fields < 10 10 G are rare, as we know 
of none. Alternatively, the dynamo model actually requires rapidly spinning systems to have 
strong magnetic fields. While core-collapse born-fast mechanisms seem to be ruled out, the 
accretion induced collapse of massive and rapidly rotating WDs to NSs might form MSPs (|22]) . 
This collapse may be able to produce the observed orbital parameters, but the large observed 
pulsar mass would require the collapsing WD to be well above the Chandrasekhar mass, and 
would also suggest that the companion should be evolved. 

Globular clusters (GCs) are known to be efficient producers of MSPs, including those in 
eccentric binaries, due to interactions between NSs and other stars or binaries in their high 
density cores. Of the ~130 known GC pulsars (|2?1) . more than 10% are in highly eccentric 
(e > 0.2) orbits. These numbers, combined with the known populations of NSs in GCs and the 
Galaxy, and the respective masses of GCs and the Galaxy, imply that GCs produce eccentric 
binary pulsars at least 1000 times more efficiently per unit mass than the Galactic disk. Further- 
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more, stellar interactions and exchanges can provide MS companions for MSPs, although those 
companions should be less massive than the most evolved MS stars currently observed in GCs 
(-0.8-0.9 Mo). 

Although GCs seem a natural formation ground for PSR J1903+0327, there is no known 
GC located near PSR J1903+0327 on the sky, nor is there any evidence for an unknown cluster 
in the 2MASS catalog (1241) . the Spitzer GLIMPSE survey data (1251) . or in our Gemini observa- 
tions. An intriguing possibility is that PSR J1903+0327 was formed in the core of a GC and 
then ejected from the cluster, possibly in the same interaction that induced the orbital eccentric- 
ity we observe today. This would have to have occurred sometime within the last — 1— 2Gyr 
(i.e. within the characteristic age of the MSP). Detailed simulations of NS interactions in GCs 
have shown that up to ~50% of recycled NSs are ejected within the ~10Gyr lifetimes of the 
GCs (26). Once the pulsar has left the cluster it drifts far away from its parent cluster on 
10 8 — 10 9 year timescales. Alternatively, the cluster could have been disrupted during orbital 
passages through the Galactic disk and bulge within the past ~1 Gyr (|27|) . Rough estimates 
based on the masses and densities of the GC system and the Galactic disk suggest (see the 
supporting online material) a 1 — 10% chance that PSR J1903+0327 originated in a GC. The 
characteristic age of the pulsar is perhaps the biggest challenge for a GC explanation as it re- 
quires that the newly-made MSP had one or more violent interactions which replaced the star 
that recycled the NS, and then was displaced enough from the GC such that the GC is no longer 
visible to us, all within ~2 Gyr. 

In a third formation scenario PSR J1903+0327 is part of a primordial hierarchical triple 
system. In such a system the inner binary evolves normally. In this case the secondary the 
secondary mass likely needs to be quite fine-tuned to to produce a highly-recycled massive 
MSP while still producing a high-mass (0.9—1.1 M Q ) WD in a wide circular orbit, although no 
such borderline systems have been observed to date. In this scenario the WD is the companion 
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seen in the timing. The third star in the system, now the MS star that we detect as the infrared 
counterpart, is in a much wider and highly inclined orbit around the inner binary. Secular 
perturbations in such a system can cause large oscillations of the inner-binary eccentricity and 
the outer-orbit inclination (so-called Kozai cycles) (l2ffl) . Especially for a cycle period resonant 
with the inner-binary relativistic periastron advance at ~ 2 x 10 6 yrs, a0.9M Q MS star ~ 120 AU 
out in a highly-inclined ~ 700 yr orbit can induce large inner-binary eccentricities (|29l) . Initial 
estimates (see the supporting online material) for the formation and survival probability of such 
a triple system suggest that a few percent of observed NS-WD binaries could be members of 
hierarchical triples. As such, it seems plausible that after finding ~50 pulsar-WD binaries, we 
have now found the first in an hierarchical stellar triple. 

A related scenario that avoids the need for formation of both a fully recycled MSP and 
a high-mass WD companion (which has never been observed before) has PSR J 1903+0327 
recycled as part of a compact inner binary in a hierarchical triple, in a configuration as recently 
suggested for 4U 2129+47 (|50l) . The MSP then ablated away its WD companion and was left 
in a 95-day eccentric orbit around the MS star we now observe. 

Further observations of PSR J1903+0327 will allow us to decide between these (or other) 
formation scenarios. A measurement of a large projected space velocity (> 100-200 km s _1 ), via 
long-term timing or Very Long Baseline Interferometry astrometry, might reflect a cluster origin 
given the high velocities of most GCs. Spectroscopic observations of the MS star will reveal 
its spectral type and metallicity, both possible indicators for or against a GC origin, and will 
show whether it exhibits the 95-day orbital motion of the pulsar. If the star is the companion, 
the radial- velocity curve will further constrain the masses of both the pulsar and the companion. 
Additionally, if the MS companion is confirmed to be more massive than ~1 M Q , it will likely 
rule out a GC origin as MS stars of that mass in clusters have already left the MS. Finally, 
long-term and higher-precision timing of the pulsar will dramatically improve the relativistic 
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parameters of the system (and therefore the derived masses) and will reveal secular changes in 
the spin and orbital parameters caused by the presence of a third star or by classical effects from 
the rotation of a MS companion. 
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Timing Parameters Assuming General Relativity 



Right Ascension (J2000) 19 h 03 m 05^79368(3) 

Declination (J2000) 03° 21' 19^2220(11) 

Spin Period (ms) 2.1499123298435(3) 

Period Derivative (s/s) 1. 879(2) x 10~ 20 

Dispersion Measure (pc cm~ 3 ) 297.537(7) 

Epoch of Period (MJD) 54280.0 

Orbital Period (days) 95.1741176(2) 

Projected Semi-Major Axis (lt-s) 105.60585(1 1) 

Eccentricity 0.436678411(12) 

Longitude of Periastron 141.65779(4)° 

Epoch of Periastron (MJD) 54063.8402308(5) 

Total System Mass, M tot (M ) 2.79(5) 

Companion Mass, M 2 (M ) 1.051(15) 

Other Parameters 

Scattering Time at 1.4 GHz (ms) 0.126(1) 

1.4-GHz Flux Density (mJy) 1.3(4) 

2-GHz Flux Density (mJy) 0.62(5) 

5-GHz Flux Density (mJy) 0.09(2) 

Spectral Index -2.1(2) 

Derived Parameters 

Galactic Longitude (J2000) 37?3363 

Galactic Latitude (J2000) -L0136 

Mass function (M ) 0.1396076(2) 

Distance (kpc) ~6.4 

Surface Dipole Magnetic Field Strength (Gauss) 2.0 x 10 8 

Characteristic Age (Gyr) 1.8 

Spin-down Luminosity (ergs s _1 ) 7.5 x 10 34 

Advance of Periastron (deg yr _1 ) 2.46(2) x 10~ 4 

Orbital inclination 78(2)° 

Pulsar Mass, M 1 (M ) 1 .74(4) 



Table 1: Measured and derived parameters for PSR J 1903 +0327. The timing parameters were 
measured using the DE405 Solar System ephemerides (37) and the "DDGR" timing model 
which assumes that general relativity fully describes the parameters of the binary system (TJ21) . 
A total of 342 pulse arrival times measured between MJDs 53990 and 54568 were fit. The 
numbers in parentheses are twice the TEMPO-reported 1-cr uncertainties in the least significant 
digit or digits quoted as determined by a bootstrap error analysis. The distance is inferred 
from the NE2001 free electron density model (14). Ensemble distance measurements using this 
model have an estimated error of 25%, although errors for individual pulsars may be larger. 
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Figure 1: Residual pulse arrival times as a function of orbital phase (mean anomaly) for 
PSR J 1903+0327 after subtraction of the best-fit timing model. The timing residuals are from 
observations made with the Arecibo telescope (blue circles), the Westerbork Synthesis Radio 
Telescope (red triangles), and the Green Bank Telescope (green crosses), and are defined as ob- 
served minus model. Top The measured timing residuals if no orbit is accounted for. The result- 
ing curve is the Roemer delay (i.e. the light-travel time across the orbit) and its non-sinusoidal 
shape shows the large eccentricity (e = 0.44) of PSR J1903+0327. The uncertainties on the 
data points would be the same as in the lower panels but the scale is different by a factor of 
10 6 . Middle The same residuals as in the top panel but with the Roemer delay and all general 
relativistic delays except for Shapiro delay from the timing solution in Table 1 removed. Bottom 
The timing residuals for the full "DDGR" timing model described in Table 1 which assumes 
that general relativity fully describes the parameters of the binary system (f32)) . The weighted 
root-mean- square timing residual shown here is 1.9 fxs. 
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Figure 2: Rotation periods, period derivatives, and orbital eccentricities (for binary pulsars) 
of pulsars in the disk of the Galaxy. The bottom face of the cube shows a plot of rotation 
period versus rotation period derivative for all Galactic pulsars. Colored points show the binary 
pulsars, projected upward from the bottom face in proportion to their orbital eccentricities. 
Square blue points are double neutron star systems, triangular green points are pulsars with 
main-sequence or massive companions, circular yellow points are pulsars with white dwarf or 
sub-dwarf companions, and the red star is PSR J1903+0327, which occupies a unique place in 
the diagram. 
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Figure 3: A Ks-band image of the PSR J1903+0327 field taken during excellent seeing con- 
ditions (0.3—0.4") with the Gemini North telescope. The red circle shows the 2-a error circle, 
with radius 0V32 (produced by the frame-tie uncertaities in right ascension and declination), for 
the position of the pulsar based on astrometric calibrations made with the 2MASS catalog. The 
star within the error circle is the possible main-sequence companion to the pulsar. 
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Materials and methods 

In the following, we describe several details of the discovery and follow-up observations of 
PSR J 1903+0327 as well as statistical arguments relating to origin in a globular cluster or a 
triple star system. 

Discovery 

The survey (5(7]) collects data during 5-min pointings for each of seven positions on the sky 
using the Arecibo L-band Feed Array (ALFA) receiver and digital spectrometers covering 100- 
MHz of bandwidth centered at 1400 MHz. The recorded data have 256 frequency channels and 
a time resolution of 64 fis. We initially detected PSR J1903+0327 in September 2006 using a 
search pipeline based on the PRESTO suite of pulsar analysis software (50 5(21). The basic steps 
of the pipeline are similar to those of the "quicklook pipeline" (5(7)), except that we process the 
data at full resolution in both time and radio frequency, we perform extensive radio-frequency 
interference excision in the time and frequency domains, and we look specifically for compact 
binary pulsars using linear "acceleration" searches (50). To remove the dispersive effects of the 
interstellar medium, the data were first dedispersed for 1272 trial dispersion measures (DMs) 
between 0— 1003 pc cm -3 . The DM trial values are spaced such that the pulse smearing due 
to interstellar dispersion is <1 ms for all DMs <600pccm 3 , thereby giving us unprecedented 
sensitivity to millisecond pulsars (MSPs) in the disk of our Galaxy. We conduct standard period- 
icity and acceleration searches by taking the Fast Fourier Transform (FFT) of each dedispersed 
time series using incoherent harmonic summing of up to 16 harmonics. In addition, single-pulse 
searches of the data were carried out which were sensitive to transient events in the time domain 
with widths in the range 64 fj,s to 0. 1 s 
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Timing 

Once a pulsar is discovered it is characterized by its position, spin parameters and any orbital 
parameters. Using the search data these parameters are poorly constrained, but they can be 
measured to much greater levels of precision by fitting a timing model of the pulsar's rotation 
to the data. Integer numbers of pulsar rotations are fitted between the times-of-arrival of the 
observed pulses; in such a "phase connected solution" every rotation of the pulsar is accounted 
for. 

Follow-up timing observations were carried out with Arecibo at ~1.5 GHz on 22 days be- 
tween September 2006 and January 2007 and at ~2.1 GHz for several days between January 
2008 and April 2008, with the Westerbork Synthesis Radio Telescope (WSRT) at ~ 1.4 GHz on 
14 days between April and December 2007, and with the Green Bank Telescope (GBT) primar- 
ily at ~2.1 GHz between December 2006 and April 2008. The Arecibo observations used the 
seven beam ALFA receiver and a Wideband Arecibo Pulsar Processor (WAPP) usually in the 
same configuration as used in the survey. The GBT observations made use of the Green Bank 
Pulsar Spigot (^O, a digital correlator which synthesized either 768 or 1536 frequency chan- 
nels covering 600 MHz of bandwidth and sampled every 81.92 /is. Individual GBT observations 
were made at 5 and 9 GHz with the Spigot in a mode with 1024 frequency channels recorded 
over 800 MHz of bandwidth and sampled every 81.92 /xs. WSRT observations used the PUMa 
II data acquisition system ($6$ to record eight 20 MHz wide dual polarization bands. Each band 
was subsequently analyzed offline by forming a 64-channel coherent filterbank leading to a final 
time resolution of 16.8 /is. 

The timing model parameters shown in Table 1 of the main article were derived from the 
best-fit model to the pulse arrival times we obtained with the above telescopes and instru- 
mentation. A timing model using a simple Keplerian orbit provides a very poor fit to the 
data (reduced-^ 2 of 29.5 for 327 degrees of freedom [DOF]). The addition of an orbital pe- 
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riastron advance (to) dramatically improves the fit (reduced-^ 2 of 2.12 for 326 DOF), giving 
<jj = 2.46(2) xl0 _4 degyr _1 , but systematic trends remain in the timing residuals. Given the 
orientation of the orbit, these remaining trends are most likely explained by the Shapiro delay: 
the delay should have an amplitude of tens of /is, an order of magnitude larger than our arrival 
time measurement precision, even though most of that signal will be covariant with the pro- 
jected semi-major axis, a sin i. Using the "DDGR" model ($7] which assumes that general 
relativity correctly describes the dynamics of the system with Keplerian orbit parameters plus 
the total system mass M to t and the mass of the companion star M 2 , the fit is further improved 
(reduced-^ 2 = 1.13 for 325 DOF) and results in no visible systematics in the residuals (Fig. 1). 
In Fig. SI, we show a x 2 contour map of the mass constraints from the DDGR model. 

In order to confirm that this timing is consistent with general relativity, the "DD" model (,5|9] 
a theory independent relativistic description of the orbital parameters, was used separately 
fitting for the orbital periastron advance, the companion mass and sine of the inclination angle 
(sini). In this fit (reduced-^ 2 = 1.12 for 324 DOF) all of the post-Keplerian parameters which 
are not specific to general relativity are consistent with the parameters from the general relativity 
specific DDGR model at the 1-2 ex level. The timing model parameters for the DD model are 
show in Table S 1 . 

Search for a pulsar companion 

Given the possibility that this is a double neutron star system, and that the companion is po- 
tentially a pulsar, we searched for pulsations from the companion using several of the early 
observations of the system taken at Arecibo. As well as searching each observation indepen- 
dently, we summed power spectra from the individual observations to increase sensitivity. We 
dedispersed each observation at the known DM of PSR J1903+0327. We then removed the 
deleterious orbital effects of the system prior to incoherently summing the power spectra from 
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FFTs. The orbital accelerations were calculated and removed by noting that the orbital period 
and eccentricity of the companion's orbit are the same as for PSR J 1903+0327, but the lon- 
gitude of periastron is different by 180°. The only unknown parameter in this calculation is 
the ratio of the semi-major axes of the pulsar and companion which is the inverse of the ratio 
of their masses. A search over a series of trial projected semi-major axes for the companion 
was therefore used, ranging from 60 to 185 light- seconds. This encompasses the range of well- 
measured neutron star masses CSlTOl). Finally the summed FFT was searched for candidates 
using the techniques described above. The search yielded no significant candidates. 

We set an upper limit on pulsed emission at 1 .4 GHz of 20 fj,Jy for a period of 2 ms and 9 fj,Jy 
for a period of 200 ms, both at a DM of 297 pc cm 3 . These upper limits correspond to 1 .4 GHz 
radio luminosities of 0.8 and 0.4 mJy kpc 2 respectively. Only two percent of all radio pulsars 
currently known have luminosities below the latter value. Our non-detection therefore excludes 
most of the observable radio pulsar population as a companion, but does not of course rule out 
a neutron star that is unfavorably beamed away from our line of sight or rotating too slowly to 
be observed as a pulsar. 

Optical and infra-red analysis 

The field of PSR J 1903+0327 was observed with the Near InfraRed Imager (NIRI) on the Gem- 
ini North telescope on July 24, 2007. Dithered series of 6 second exposures were obtained in the 
J, H and K$ filters, amounting to a total exposure time of 10 minutes in each filter. All images 
were corrected for dark current and flatfielded using skyflats. The non-uniform sky distribution 
was removed by substracting a sky frame constructed of the unregistered science images to re- 
move the contributions of stars. After these corrections, the images taken with the same filter 
were registered and averaged. 

Astrometry was done relative to the 2MASS catalogue (Sll), as the small field of view of 
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NIRI (2' x 2') contains no astrometric standards from the UCAC2 catalogue (SI 2). A total of 71 
2MASS stars coincided with the registered and average _K"s-band image and 51 of these were 
not saturated and appeared steller and unblended. After iteratively removing 9 outliers, the final 
astrometric calibration has rms residuals of 0'/086 in right ascension and 0'.'094 in declination. 
A comparison of the positions of 281 UCAC2 standards coinciding with 2MASS stars within 
10' from the position of PSR J1903+0327 shows no significant shift between the positions in 
both catalogues. 

Taking the rms uncertainty on both coordinates as a measure for the l-a uncertainty in the 
astrometric calibration, we find a single star inside the 0'.'13 error circle on the pulsar position at 
aj2ooo = 19 h 03 m 05!796(6) and 5 J2 ooo = 03°27 / 19'.'19(9), where the uncertainty on the position 
is the quadratic sum of the uncertainty in the astrometry and the intrinsic positional uncertainty 
of the star in the image (0'.'02 in each coordinate). This position if offset from the pulsar position 
by 0"034(88) in right ascension and 0"032(96) in declination. 

Photometry was also done relative to the 2MASS catalogue. Magnitude offsets were deter- 
mined and outliers were iteratively removed, leaving 30 to 40 stars for the photometric calibra- 
tion. The calibrations showed no significant dependence with star color. We find that the object 
in the error circle has J = 19.22(9), H = 18.41(10) and K s = 18.03(9). 

To estimate the reddening towards PSR J 1903+0327, we used red clump stars to trace the 
reddening as a function of distance ( Sd3\i . Again using the 2MASS catalogue, we selected 1400 
stars within 5' from PSR J 1903+0327 and determined the J — Ks colour of the red clump stars at 
different K s magnitudes using the formalism of ( SlT4\i . At the DM-distance of PSR J1903+0327, 
d = 6.4(16) kpc (assuming a 25% uncertainty), we constrain the reddening at Ay = 4.9(11). 
For this distance and reddening, the intrinsic magnitudes of the star inside the error circle are 
Mj = 3.8(11), M H = 3.5(10) and M Ks = 3.5(10). 

Though the metallicity and age of the star are unknown, stars of solar metallicity, ages of 
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1 Gyr or less and masses between approximately 0.8 and 1 .3 M Q will have similar absolute mag- 
nitudes (S15). Stars older than 1 Gyr will have evolved and will match the absolute magnitudes 
for lower masses. At 10 Gyr and solar metallicity the mass range decreases to approximately 
0.75 to 1.1 M . For a mass of 1.05 M as determined from the radio timing, a solar metallicity 
star will match the observed absolute magnitudes if the age is less than 6 Gyr. 

At the apparent magnitudes, the star is bright enough for optical or infrared spectroscopy to 
determine radial velocity variations which would unambiguously confirm the star as the binary 
companion to PSR J1903+0327. 

Globular cluster calculations 

While PSR J1903+0327 is unprecedented in the Galactic plane several MSPs in eccentric bi- 
naries have been seen in globular clusters (GCs). GCs are known to be excellent breeding 
grounds for MSPs, with formation rates per unit mass ~100 times that in the disk CSfiol §Th . 
PSR J1903+0327 may have formed in a GC which was later disrupted, or escaped from an 
existing GC ~ 10 9 yr ago such that they are no longer in the same part of the sky. Many GCs 
may have been completely disrupted via gravitational interactions with the disk and bulge over 
the course of many orbits in the Galactic potential. As a result, their stars become part of the 
Galactic spheroid population (an approximately spherical distribution of older stars distributed 
in an extension of the central Galactic bulge, with a diameter of ~10kpc). Estimates suggest 
that over half the spheroid mass could have come from such disrupted GCs CSfTffi). Such dis- 
ruption should require several Galactic orbits of the GC, the periods of which are ~0.5— 1.0 
Gyr. However, in the particular case of PSR J 1903+0327, we are constrained by the relatively 
small characteristic age of the system. Since a characteristic age is usually an upper limit on the 
true age of a radio pulsar (3?]), any model involving GC disruption requires a recycled pulsar 
to exchange its companion and acquire an eccentric orbit and its parent GC to then be totally 
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disrupted in less than ~2 Gyr. This is challenging given the orbital periods of GCs. 

By comparing the densities of the spheroid and disk we can estimate the probability of 
PSR J 1903+0327 being found in the disk but actually being associated with the spheroid pop- 
ulation. Using estimates of the stellar densities at low Galactic latitudes of both the spheroid 
(^O and disk populations, we find p S p h c r oid/Pdisk = 4.7xl(T 5 M Q pc~ 3 /4.3x 1(T 2 M Q pc~ 3 ~ 
1 x 1CT 3 . Assuming that half of the spheroid mass originated in disrupted GCs, the <2 Gyr age 
of PSR J1903+0327 implies that <20% of that mass could have come from GCs disrupted since 
the recycling of the pulsar. Furthermore, if we assume that those GCs produced highly eccentric 
binary MSPs ~10 3 times more efficiently per unit mass than the disk, we can crudely estimate 
that the probability that PSR J1903+0327 came from a disrupted GC is <10%. 

We can make a similar estimate of the probability that PSR J 1903 +0327 was ejected from 
a GC. The number of GC-ejected PSR J1903+0327-like systems should be roughly the ratio 
of the GC system mass (10 7 M Q ) to the disk mass (10 11 M ) times the 10 3 formation efficiency 
factor for highly eccentric binary MSPs times the ~50% fraction of pulsars which are eventually 
ejected from GCs ($21\ . The resulting ~5% probability again indicates that some GC formation 
mechanism for PSR J1903+0327 is plausible. 

Hierarchical triple system calculations 

Hierarchical stellar triple systems are common, and such a system may explain the eccentricity 
we observe in PSR J1903+0327. To date no pulsars have been found in other stellar triples. 
Hierarchical triples contain an inner binary (in this case the MSP and the ~1.05 M Q companion 
seen in timing at the 95-day orbit) with a third object orbiting at greater distance (in this case, 
the main sequence star observed in the infra-red). Here we assess the plausibility of forming 
such an end state. 

Constraints on the formation of the inner binary are: (i) to form a 2.15-ms MSP there must 
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have been a long period of stable mass transfer to spin up the neutron star; (ii) to produce the 
wide 95-day orbit the white dwarf progenitor cannot be too massive, to avoid common-envelope 
(CE) evolution and spiral-in; (iii) the star must have been massive enough to leave a ~ 1.05 M Q 
white dwarf after final mass transfer ($22\. We therefore consider an initial system containing 
a central binary formed by a 8 M Q primary; a secondary that is just below CE mass, we shall 
here assume 3 M ; and a 0.9 M third star. 

A survey of the multiplicity of massive stars in the Orion Nebula ($23\i suggests that 80% 
of spectroscopic binaries with OB-star primaries have a third visual companion. These systems 
must survive the supernova that creates the pulsar: the inner binary will remain bound given 
suitable mass-loss and kick parameters (S24). The outer companion can remain bound CSl25l) 
if its orbital velocity is comparable to the new system velocity of the inner binary (a few tens 
of kms -1 , ($26\i). The multiplicity survey ($23\i finds ~5-10% of third companions are close- 
in enough to remain bound in the supernova. The semi-major axis of this outer orbit may 
then increase as mass is lost from the system during the pulsar spin-up phase ($2~7\i. In this 
hierarchical triple the outer companion cannot unbind the inner binary, but its orbit evolves 
towards a Kozai-cycle resonance with the inner-binary relativistic periastron advance period 
(S28). Our above rough estimates of the various fractions combine to a roughly 4% chance that 
a given evolved NS-WD binary is part of a hierarchical triple. As there are currently ~50 NS- 
WD systems known in Galactic-disk binaries, this simple estimate is consistent with the idea 
that PSR J1903+0327 is the first in such a triple. We do note that this formation scenario still 
faces the difficulty that thus far no other such systems with both a highly recycled pulsar and a 
wide-orbit very massive WD companion have been observed. 

Two other probability arguments for the detection of this system are not directly related to 
the formation plausibility. First, for a system with masses 1.74, 1.05 and 0.9 M Q , a relatively 
high inclination of the third star's orbit i > 45°and some eccentricity in the outer orbit, the 
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inner-binary eccentricity is higher than e = 0.44 for about 20% of the oscillation time (S28). 
Second, for the line-of-sight acceleration from the third star to be less than our measured P this 
third star must currently be close to plane of the sky, the a priori probability of which is ~5%. 
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Timing Parameters Using the DD Model 



Right Ascension (J2000) 19 h 03 m 05!79368(3) 

Declination (J2000) 03° 27' 19"2220(1 1) 

Spin Period (ms) 2.1499123298435(3) 

Period Derivative (s/s) 1. 879(2) x 1(T 20 

Dispersion Measure (pc cm -3 ) . . 297.537(7) 

Epoch of Period (MJD) 54280.0 

Orbital Period (days) 95.1741 176(2) 

Projected Semi-Major Axis (lt-s) 105.593459(8) 

Eccentricity 0.43667838(5) 

Longitude of Periastron 141.652477(3)° 

Epoch of Periastron (MJD) 54063.8402310(6) 

Advance of Periastron (deg yr _1 ) 2.46(3) x 10~ 4 

Sine of Orbital Inclination 0.966(10) 

Companion Mass, M 2 (M ) 1-3(2) 



Table SI: Measured and derived parameters for PSR J1903+0327. The timing parameters 
were measured using the DE405 Solar System ephemerides (S29) and the "DD" timing model 
which uses a theory independent relativisitc model to describe the parameters of the binary 
system (5[P] A total of 342 pulse arrival times measured between MJDs 53990 and 54568 
were fit. The numbers in parentheses are twice the TEMPO-reported l-a uncertainties in the 
least significant digit or digits quoted as determined by a bootstrap error analysis. 
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M T0T (Solar Masses) 



Figure S 1 : A x 2 -map showing the 1 and 2-a confidence regions for the total system mass (M tot ) 
and companion mass (M 2 ) based on the timing solution presented in Table 1 after doubling the 
nominal TEMPO errors. Those errors were within 20% of the bootstrap error estimates. The 
grey region in the lower right portion of the figure is excluded by the mass function. The grey 
lines show inclinations i of 70, 75, 80, and 85°. For PSR J1903+0327, M tot is best constrained 
from the measurement of the relativistic advance of periastron uo while M 2 is constrained by a 
measurement of sin % via the Shapiro delay (see Fig. 1). 



28 



